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We report experimental results for a boundary-mediated wavenumber-adjustment mechanism and 
for a boundary-limited wavenumber-band of Taylor- vortex flow (TVF). The system consists of fluid 
contained between two concentric cylinders with the inner one rotating at an angular frequency 
Q. As observed previously, the Eckhaus instability (a bulk instability) is observed and limits the 
stable wavenumber band when the system is terminated axially by two rigid, non-rotating plates. 
The band width is then of order e 1 / 2 at small e (e = f2/f2 c — 1) and agrees well with calculations 
based on the equations of motion over a wide e-range. When the cylinder axis is vertical and 
the upper liquid surface is free (i.e. an air- liquid interface), vortices can be generated or expelled 
at the free surface because there the phase of the structure is only weakly pinned. The band of 
wavenumbers over which Taylor- vortex flow exists is then more narrow than the stable band limited 
by the Eckhaus instability. At small e the boundary-mediated band-width is linear in e. These 
results are qualitatively consistent with theoretical predictions, but to our knowledge a quantitative 
calculation for TVF with a free surface does not exist. 



PACS Numbers: 47.54.+r, 47.20.-k, 47.32.-y 



I. INTRODUCTION 



Many physical systems will undergo a transition from 
a spatially-uniform state to a state with spatial variation 
when they are driven sufficiently far from equilibrium by 
an external stress. This spatial variation is known as a 
" pattern" . Typical of such systems are Rayleigh-Benard 
convection (RBC) jl| and Taylor- vortex flow (TVF) @. 
The present paper will focus on TVF, which is a flow pat- 
tern of a fluid confined between two concentric cylinders. 
When rotating the inner cylinder about its axis, the ini- 
tial instability which occurs as the angular frequency O 
is increased is a transition at f2 c from a uniform (circu- 
lar Couette) state to the TVF state. TVF is periodic in 
the axial direction and has full rotational symmetry in 
the azimuthal direction. In the axial direction the pat- 
tern is characterized by a wavenumber k — 2it/\, where 
A is the wavelength. One wavelength consists of a vor- 
tex pair since adjacent vortices have opposing circula- 
tion. For fi > O c , TVF solutions to the equations of 
motion exist over a band of wavenumbers which depends 
on e = f2/f2 c — 1 and on the boundary conditions in the 
axial direction. 

We refer to patterns with non-trivial variation only 
in one direction (such as TVF) as one-dimensional pat- 
terns. It is well known that the band of wavenumbers of 
one-dimensional patterns of finite length, with the pat- 
tern phase pinned at each end, is limited by the marginal 
curve where perturbations of the uniform state first ac- 
quire a positive growth rate. In TVF with rigid, non- 



rotating ends a large-amplitude Ekman vortex forms ad- 
jacent to each end and is associated with phase pin- 
ning. However, the resulting solutions of such a system 
are not all stable. The (more narrow) stable band of 
states is limited by the long- wavelength Eckhaus instabil- 
ity. §,||g,@,|,|^Jl|||,||,|| The Eckhaus instability 
is a bulk instability which manifests itself in the system 
interior where one pattern wavelength (one pair of Taylor 
vortices) is either gained or lost, depending on whether k 
is larger or smaller than the stable band limited by fcs(e). 
One of the successes in the study of pattern-forming sys- 
tems is the agreement between the calculated and mea- 
sured fcg(e) for TVF which has been found with three 
different ratios of the cylinder radii. Pj7|,p|, ^0| , pT[ 

Theoretically, phase pinning (and thus the Eckhaus 
instability) at small e is expected when the boundary 
conditions at the system ends z = 0, L correspond to a 
large amplitude A(Q) = A(L) = Ao of the velocity field, 
say Aq — 0(1), while in the system interior the ampli- 
tude A(z) is small, say 0(e 1 ^ 2 ). This situation closely 
corresponds to the TVF system with rigid ends where 
the influence of the Ekman vortex can be approximated 
by this boundary condition. |i^|J^,[lo|] The Eckhaus-stable 
band then has the same width as the band of stable states 
for the infinite system, but the number of states is 
finite since only discrete states with an integer number 
of vortices N = 2L/X = Lk/ir can occur. 

It was shown theoretically by Cross et al. JL7 18 
(CDHS) and discussed from various viewpoints by oth- 
ers @|0[|l]j22|j2||,|,|ll) that a reduction of A to A a = 
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Ae 1 / 2 with A < 0(1) leads to a qualitative change of the 
width of the wavenumber band over which solutions ex- 
ist. Solutions corresponding to such boundary conditions 
are known as type I solutions. The limit of their ex- 
istence is determined by phase slip at the boundary. The 
solutions appear to be stable over their entire existence 
range. The band limits kt and k^ at large and small k 
respectively vary linearly with e when e is small, yielding 
kfr = A~e and = X + e where A are constants which 
depend on the particular system. Thus a dramatic reduc- 
tion of the band at small e from 0(e 1 ^ 2 ) (in the Eckhaus 
case) to 0(e) is predicted. Solutions corresponding to 
A > 0(1) are called type II solutions; Q they exist over 
a wide range inside the marginal-stability curve and are 
stable over the entire Eckhaus-stable band. 

The reduction of the band and the occurrence of phase 
slip near the sidewall was verified by experiments on RBC 
|^J,^5|, on the buckling of plates subjected to a compres- 
sional stress |^,^|, and on Benard-Marangoni convec- 
tion p8[ . However, this early work was mostly qualita- 
tive and for relatively large e. Qualitative confirmation 
of a boundary-reduced wavenumber band was obtained 
also from numerical calculations for narrow RBC systems 
of finite length. ||,|o) 

More recently, Mao et al. investigated boundary- 
limitation of the wavenumber-band in a long, narrow, 
and very thin (two-dimensional) film of a smectic-A liq- 
uid crystal which undergoes convection when driven by 
an electric field. |31| The convection rolls form a one- 
dimensional pattern. Measurements showed that the flow 
velocity of the rolls is strongly reduced near the ends of 
the system, suggesting that the boundary conditions cor- 
respond to a value of A much less than one. Thus one 
would expect the boundary-mediated mechanism to come 
into play. Indeed it was observed in the experiment that 
convection rolls were always lost or gained at the system 
ends and never in the interior. For this system the lo- 
cation of the Eckhaus stability limit was estimated from 
the location of the neutral curve, which in turn was cal- 



culated from the equations of motion |32|. Thus a com- 
parison of the reduced boundary-limited band with the 
approximate location of the wider Eckhaus band of the 
infinite system is possible. The results include data over 
the range < e < 1. They convinvingly demonstrate 
the reduced band width and are consistent with a linear 
dependence of k^ — k c on e in the range e < 1 {k c is the 
critical wavenumber which forms for e = 0). 

Here we report experimental results for the wavenum- 
ber band of TVF. Using rigid, non-rotating ends, we re- 
produced the well known Eckhaus mechanism and sta- 
bility range. [0J|]. We then investigated a system with a 
free liquid-air interface which terminated the fluid axially 
in an apparatus oriented vertically with its axis. Since 
the influence of the viscosity of air and surface-tension 
effects presumably were small, one would expect such a 
surface to be approximated reasonably well by free-slip 
boundary conditions in the theory, and these would lead 



to an amplitude A$ equal to that in the bulk (i.e. A ~ 1). 
Semi-quantitative verification of this was obtained from 
the flow visualizations to be presented below. Thus, we 
expect this TVF case to be on the borderline where it 
would be difficult to say a priori whether the observ- 
able band is limited by the Eckhaus (bulk) instability 
or whether the existing band is limited by phase slip at 
the boundary. We demonstrated experimentally that the 
free liquid-air interface provides sufficiently weak phase 
pinning to permit phase slip. Thus vortices were gained 
or lost at that surface rather than in the bulk. This 
mechanism led to a narrower band than the Eckhaus- 
stable band, with boundaries at fc b (e) which at small e 
were linear in e as expected theoretically. However, these 
boundaries did not pass through k c and e = 0. Instead 
there remained a small gap. It is not clear whether this 
gap was due to the finite length of our system Jll| or 
whether it was associated with the marginal boundary 
conditions (A ~ 1). 

Although it is peripheral to the main purpose of this 
paper, we note that the radial velocity component at the 
free surface was almost always found to be in the out- 
ward direction. This is consistent with free-slip bound- 
ary conditions, where the centrifugal force is unopposed. 
It differs from the rigid boundary near which normally 
the flow is inward because the rigid boundary condition 
of vanishing velocity can not sustain the radial pressure 
gradient required by outward flow. Thus, with one rigid 
and one free boundary, our system normally contained an 
odd number N of vortices. However, anomalous states 
with inflow at the free surface (and thus an even N) were 
also encountered on rare occasions. 

In Sect. H of this paper we describe the experimental 
apparatus and procedures, and the methods of analysis 
which were used. The experimental results for the insta- 
bility mechanism and the stability band are presented in 
Sect. III. A brief discussion of the results and some ideas 



for future work follow in Sect. IV 



II. APPARATUS AND PROCEDURES 

A schematic diagram of the apparatus is shown in 
Fig. [|. The Couette- Taylor column was similar to the 
one described by Dominguez-Lerma et al. 0. It con- 
sisted of two concentric, straight cylinders. The outer 
one was made of plexiglass. The inner cylinder had 
an aluminum core which was clad with delrin. The 
radii were r\ = 18.68 and r-i = 25.44 mm, yielding 
a radius ratio r\ = r\jr% = 0.734 and a radial gap 
d = T2 — r\ = 6.76 mm. A computer-controlled step- 
per motor rotated the inner cylinder with an angular fre- 
quency Q — 27r/, and the outer cylinder was at rest. 
The system was temperature controlled to about ± 20 
mK by water from a Neslab refrigerated circulator. This 
water flowed through a transparent jacket surrounding 
the outer cylinder. 
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FIG. 1. Schematic diagram of the apparatus. 



The gap between the cylinders contained two non- 
rotating delrin collars, one located near each end of the 
apparatus. The collars had a flat surface facing the fluid 
and orthogonal to the cylinder axis. The axial position of 
the upper collar was adjustable by means of two 2.5 mm 
diameter stainless-steel rods which passed through seals 
in the upper end cap of the apparatus and attached to 
the collar. We were able to create either a rigid or a free 
(liquid/air) surface at the top of the column by adjusting 
the liquid level and/or the top collar position. This also 
allowed for changes of the aspect ratio L = H/d, where 
H is the distance in mm between the rigid boundaries 
or between the bottom boundary and the free surface. 
Typically, we used L ~ 40 to 50. 

The fluid level in the Taylor column could be adjusted 
by changing the vertical position of the reservoir shown in 
Fig. This reservoir was attached to a vertical transla- 
tion stage which could be moved with a second computer- 
controlled stepper motor. 

For the working fluid we used a mixture of 50% glyc- 
erol, 48% water, and 2% Kalliroscope for flow visualiza- 
tion. This mixture had a kinematic viscosity of about 
0.07 cm 2 /s, yielding a gap-diffusion time t u = d 2 jv ~ 7 
sec. The axial diffusion time Y?t v then is about 4 hours 
for the typical values L ~ 45 used by us. Typical ax- 
ial equilibration times required to reach an axially uni- 
form state are determined by phase diffusion and depend 
on the wavenumber. For small e and k near k c they 
are expected f^U^l to be about 0.2L 2 t ll ~ 50 minutes. 
Near the Eckhaus boundary they are expected to be- 
come much longer because the phase diffusivity vanishes 
at that boundary. We depended upon the measured time 
evolution of the phase of our system to decide whether 
the data were quasi-static. Amplitude equilibration after 
a change of / is generally much faster and typically oc- 
curs on a time scale of Tot„/e with |35| to = 0.076, which 
for e > 10~ 2 is only about a minute. 
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FIG. 2. Illustration of the image analysis. The top is an ex- 
ample of a single image (a "snapshot"). The middle is a plot 
of the image amplitude, averaged in the narrow (azimuthal) 
direction, as a function of the axial direction in units of cam- 
era pixels. The bottom is the Fourier transform of the am- 
plitude (after length was scaled by the gap d) and shows the 
wavenumbers of the Taylor vortices. 



A CCD camera was used to digitally record an im- 
age of the Taylor-vortex pattern. Typically, we took 
images with a width (in the azimuthal direction) of 30 
pixels and a length (in the axial direction) of 534 pix- 
els. An example is shown in Fig. |^a. Division by 
a reference image taken below the onset of TVF was 
used to reduce inhomogeneities due to uneven illumina- 
tion. The pixels were averaged over the width to pro- 
duce a one-dimensional record of the TVF state as a 
function of the axial position. Temporal sequences of 
such one-dimensional records were used to prepare a two- 
dimensional space-time "image" which showed the evo- 
lution of the TVF state as a function of time and/or 
inner-cylinder frequency (see for instance Figs. ^, ^, and 
^ below). In order to avoid effects associated with the 
cylinder ends (such as the strong Ekman vortex at the 
rigid bottom boundary), only a central segment of length 
256 pixels of each one-dimensional record was used for 
further analysis. The section used is shown in Fig. |]a by 
the black bar. An example is plotted in Fig. |^b. The vari- 
ation of the signal amplitude with axial position is due 
primarily to effects associated with uneven illumination 
which remained even after image division. 
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FIG. 3. An illustration of the determination of tt c . The top 
image is a space-frequency (or time) plot with an expanded 
frequency scale near onset. The middle gives the signal power 
determined from the Fourier transform (see Fig. ^|). The bot- 
tom shows the corresponding wavevector determinations. 

The data were multiplied by a Harming window func- 
tion and then Fourier transformed. The square of the 
modulus of such a transform is shown in Fig. |p as a 
function of the wavenumber fc. It is apparent that the 
second harmonic is significantly stronger than the fun- 
damental. This is expected from the Kalliroscope flow 
visualization since outflow and inflow boundaries of vor- 
tices have a similar appearance. We computed the first 
moment of the second harmonic, using typically only 
the three points nearest the maximum, to determine the 
wavenumber of the TVF structure. To determine the on- 
set frequency f c , we also determined the power under the 
second-harmonic peak. 

In order to determine f c , a TVF state was usually pre- 
pared first by slowly increasing the cylinder frequency / 
to a value above the critical value f c , and by waiting for 



equilibration of this state. Then / was stepped down in 
increments of 2 mHz, with an equilibration time of 100 
to 600 seconds after each step. On average this corre- 
sponds to a ramp rate (3 = {t v j f c )(df /dt) < 10~ 4 (here 
t is in sec). The image in the top of Fig. |^ shows the 
evolution of the axially- varying amplitude for the section 
of such a run very near f c as a function of time and thus 
of frequency. Here the individual time/frequency steps 
are resolved. For this case there was a rigid boundary 
at each end. The Ekman vortices and their influence on 
the internal TVF structure are apparent. However, the 
system is sufficiently long that the ends do not signif- 
icantly influence the determination of f c in the interior. 
The middle part of Fig. || gives the power of the second 
harmonic of the Fourier transform (determined as dis- 
cussed above) as a function of / over the same range as 
the image above it. For the example in the figure it is 
easy to see that f c = 646 ± 2 mHz. Increasing the fre- 
quency from below f c gave the same results within the 
quoted uncertainties. The bottom part of Fig. || shows 
the wavenumber derived from the same run. One can see 
that fc is determined unambiguously by the data at all 
frequencies greater than f c . 
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FIG. 4. Illustration of the experimental protocol (see text 
for details) . The solid line is a smooth curve passing through 
our results for the boundary-mediated stability limit. The 
dash-dotted curve is the Eckhaus stability limit. 

Two examples of an experimental protocol designed to 
determine the boundary-selected stability limits are il- 
lustrated in Fig. [|. When the large- (small-) k stability 
boundary was to be investigated, a liquid column with 
a free surface and with a large (small) L = Lq was first 
prepared. The examples in Fig. |f| correspond to Lq =61 
(37) . An initial state with k near k c was then prepared by 
slowly ramping / from below to above f c . A state with 
N = 61 (37) vortices resulted. Next the frequency was in- 
creased relatively quickly to two to four times f c (points 
Ai and A2 in the figure), generally without changing the 
number of vortices. The wavenumber of the system was 
then adjusted to a desired starting value [point Bi (B2) 
in the figure] by lowering (raising) the liquid level in the 
TVF column. This was accomplished by decreasing (in- 
creasing) the vertical reservoir position. When the reser- 
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voir position was changed sufficiently slowly, no loss or 
gain of vortices occurred provided the structure remained 
well within the stable wavenumber band. Typically we 
changed the reservoir position at the rate dz/dt ~ 10~ 3 
cm/s, corresponding to (t v / d)(dz / dt) ~ 1CT 2 . Typically 
one or two hours were required to reach points Bi or B2. 
In the examples of Fig. || B± (B2) corresponded to L = 
45.9 (51.3). 

Having established states Bi or B2, the frequency was 
reduced in steps of 2 mHz at time intervals of typically 
100 to 600 seconds, corresponding to effective ramp rates 
near j3 = — 10~ 4 . Before each new step, an image was 
taken. The circles in Fig. | illustrate idealized sequences 
of instability points which are encountered in this pro- 
cess. Actual experimental runs will be illustrated in the 
next Section by space-time images and by plots of the 
wavenumber as a function of time or frequency. 

The procedure for determining the Eckhaus boundary 
was similar, except that at the state indicated by Hi and 
B2 in Fig. [|the rigid collar bounded the fluid at the top. 
Sudden-start protocols were used as well at times to 
quickly prepare initial states with k different from k c . 
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FIG. 5. Illustration of the two mechanisms for wavenumber 
limitation. 



III. RESULTS 



Figure || gives space-time images which illustrate the 
two different mechanisms. The left part is a run with 
rigid top and bottom ends, and with k < k c . The Ek- 
man vortices and their influence on the amplitude of the 
Taylor vortices in the interior but near the two ends are 
apparent. The instability mechanism led to the gain of 
vortex pairs in the system interior. Three transitions 
are clearly visible, although the third was almost im- 
mediately followed in the sample interior by the axially- 
uniform Couette state. 

A similar experiment with a free surface yielded a very 
different result. This is illustrated in the right part of 
the figure. In this case the original state (near the right 
of the image) has k > k c . Five or more transitions can 
be seen, but vortex pairs were lost at the free top surface 
rather than in the interior. 

The phase slip at the free surface is illustrated in 
greater detail in Fig. ^. The left example is for k < k c . 
In this case vortices were added at the free surface, so as 
to increase k. The right one is for k > k c where vortices 
were expelled so as to decrease k. Figure ^| also illustrates 
the nature of the boundary condition at the free surface. 
The TVF amplitude is essentially constant as a function 
of z, with at most a very small enhancement (depression) 
visible for k < k c (k > k c ) as the Couette state is entered. 
The very small difference, if any, between the two cases 
shown probably is due to varying conditions with height 
of the cylinder surfaces rather than a dependence on k. 




FIG. 6. Detailed view of the losses (right image) and gains 
(left image) of vortices at the free upper surface. 

Figure [7|a gives the wavenumber as a function of / for a 
run which started with k < k c . The parameters are those 
used in Fig. || to illustrate the small-A: experiments. Here 
L = 51.3, and initially (at the right of Fig. 0a) there were 
N = 37 vortices. A gain of a vortex pair corresponds to 
SN = 2 and should yield Sk = {n/L)SN = 0.122 if the 
anomalous width of the Ekman vortex at the bottom 
boundary is neglected. The average measured step size 
was also 5k = 0.122, consistent with the expected effect 
of a pair gain. It is clear from Fig. fj]a that most of the 
transitions are pair gains, thus maintaining an odd num- 
ber of vortices in the system. This is consistent with the 
preference for outflow at the free and inflow at the rigid 
surface. However, there are occasional exceptions. One 
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of them can be seen in Fig. [?]a near / = 722 mHz, where 
two successive single-vortex-gains (TV = 47 to 48 to 49) 
resulted in smaller values 8k ~ 0.058. 
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FIG. 7. The wavenumber k as a function of the in- 
ner-cylinder frequency / . Temporally, the experiment pro- 
ceeded from right (large /) to left (small /). (a): k < k c . (b): 
k ^> kc ■ 

It can be seen that the system equilibrates surprisingly 
fast after each transition. This is qualitatively consistent 
with calculations of the phase diffusivity D(k,e), which 
show a strong asymmetry as a function of k at constant 
e with a maximum at relatively small k < k c . j3?J Equi- 
libration is indeed much slower for k > k Cl as can be 
seen in Fig. 0b. For this case L = 45.9 and N — 61 at 
the beginning of the run. After each transition there is a 
noticeable slow relaxation of the average wavenumber in 
the section of the system used for the analysis. Nonethe- 
less, reasonable estimates of the stability boundary can 
be obtained from the experiment. In this case the steps 
have an average value 8k = 0.141. For pair losses and 
this L one would expect 0.136, in good agreement with 
the experiment. 

Figure ^ summarizes all our determinations of the sta- 
bility boundary with a free surface as solid circles. They 
were obtained from numerous runs with different values 
of L over a time period of about 4 months, each run giv- 
ing a few points on one of the boundaries. They scatter 
somewhat more than we anticipated on the basis of our 
resolution for e and kf. We believe that there may be a 
contribution to the scatter from irreproducibilities from 
run to run in the nature of the air-liquid interface. This 
interface could for instance be influenced by contamina- 



tions of the cylinder surfaces by the Kalliroscope, and this 
contamination could vary with L; but we do not have any 
direct information about this. For reference, we show the 
theoretical Eckhaus boundary as a solid line and the ex- 
perimental determinations of the Eckhaus boundary by 
Dominguez-Lerma et al. as open circles. Regardless of 
their scatter, it can be seen that the free-surface-limited 
data lie well inside the Eckhaus-stable band over the en- 
tire e-range of the experiment. 
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FIG. 8. Results for the wavenumber-band limits in the 
presence of a free boundary (solid circles) . Also shown (open 
circles) are the results from Ref. Jt) for the Eckhaus boundary 
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FIG. 9. Results at small e for the wavenumber-band limits 
in the presence of a free boundary (solid circles) . Also shown 
(open circles) are the results from Ref. Q for the Eckhaus 
boundary and the theoretical Eckhaus boundary (solid line) 
from Ref. §. 

At small e the data are consistent with a linear de- 
pendence on e, as shown by the dashed lines in the fig- 
ure. The small-e portion of the results is shown once 
more on expanded scales in Fig. ^|. The linear depen- 
dence of k s upon e is clearly seen. The slopes k c de/dk^ 
of the boundaries are within error symmetric and equal 
to about ±1.7. However, the straight lines through the 
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data do not pass through the origin. Instead, there is 
a small gap near k c which reaches approximately from 
Sk-/k c = -0.01 to 5k+/k c = 0.03. Whereas the width 
8k + /k c — 5k~ /k c ~ 0.04 of the gap has an uncertainty of 
only about 0.01, the values of 5k + /k c and 8k~ /k c them- 
selves are less well known because of possible small sys- 
tematic errors in k. Thus a gap symmetric about k c is 
not ruled out. 



IV. DISCUSSION 

We presented experimental results which demonstrate 
that TVF terminated at one end by a liquid-air inter- 
face has a band of solutions which is more narrow than 
the Eckhaus-stable band of the system with rigid non- 
rotating ends. The band width is limited by phase slip 
at the free surface. At small e, the band limits fc^(e) vary 
linearly with e. This is in agreement with general theo- 
retical predictions, but to our knowledge specific calcula- 
tions for TVF have not been carried out. The boundaries 
fc^(e) do not pass through k c at e = and instead leave 
a small wavenumber gap Sk + — 5k~ at the neutral curve 
which was not theoretically predicted. 

At this point we can only speculate about the origin 
of the gap. Perhaps it is asociated with the particular 
boundary conditions which pertain to our physical sys- 
tem. As discussed above, we expect TVF with a free sur- 
face to be at the borderline between systems for which 
the solution band is limited by phase slip at the bound- 
ary (type I solutions) and those which have stable so- 
lutions over the entire bulk Eckhaus instability (type II 
solutions). Perhaps the transition from type I to type II 
solutions takes place by a gap opening near k c , although 
a more likely course of events would involve going to 
zero (i.e. a vanishing of the slopes of the straight lines in 
Fig. |). 

Perhaps a more likely explanation can be found in the 
finite length of the system. For the Eckhaus boundary 
of the finite system it is known that a gap opens up [^||^] 
near k c which, for periodic boundary conditions, |^] has 
the width ±Sk/k c = ±ir/Lk c . Thus, any state with 
—w/Lk c < (k — k c )/k c < ir/Lk c is stable all the way 
down to the neutral curve where its amplitude vanishes. 
For L ~ 50 this gap in the Eckhaus boundary is shown 
in Fig. |9| by the two short dashed lines. One sees that 
the gap Sk^/kc derived from our data for kr(e) is of the 
same size; but as yet there is no theoretical prediction 
for finite-size effects on k^(e). 

Finally we comment on possible future experiments. 
Clearly the small-e, small- (fc — k c ) range needs to be ex- 
amined by more detailed experiments to gain further in- 
sight into the nature of the gap discussed above. In part 
this can be done by varying the system length to see 
whether the gap is proportional to L^ 1 . Another in- 
teresting possibility for future work will be an attempt 
to alter the boundary conditions, thus effectively chang- 



ing A. It may be possible to do this by placing a less 
dense but perhaps more viscous fluid above the glycerol- 
water solution. Since the flow field of the Couette state 
is independent of the viscosity, one would not expect an 
Ekman vortex to be generated. However, the larger vis- 
cosity would lead to subcritical conditions in the added 
fluid when Q = Q c in the working fluid. The ensuing 
Couette state, located physically above the TVF state, 
should thenjeduce the TVF amplitude. This would cor- 
respond to A < 1 and should yield a more narrow range 
of solutions (presumably larger |A ± |) than those found in 
the present work. 
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